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The Thrust Optimised Parabolic nozzle 


Introduction 

In the early 1950’s, rocketeers attempted to devise ‘the perfect nozzle’: one that would cause 
the least thrust losses. They used the semi-algebraic Method of Characteristics (of supersonic 
flow) to devise such a nozzle, for whatever nozzle expansion ratio € was required. 

Unfortunately, at large expansion ratios, this nozzle was far too long, resulting in excess 
nozzle mass. 

Rao in America and Shmyglevsky in Russia found a way to modify this method in order to 
produce an optimum nozzle that was much shorter: it resembled a church bell and was hence 
known as a ‘bell nozzle’. 


For ease of use, Rao measured the lengths of his resulting nozzles as fractions of the length 
of a standard 15 degree half-angle conical nozzle which has length: 

j (Ve-i) R t 

L N cone = - -—— equ. 1 where R t is the radius of the throat. 

“ tan^ibj 


So an ‘80% bell’ would have a length of 0.8 times this length. 


Rao was clearly also schooled in traditional geometry. He was able to spot that the bell part of 
his bell nozzles could be approximated very closely by a skewed parabola allowing us to 
quickly sketch his nozzles with negligible loss of thrust performance. 

These sketched approximations are known as Thrust Optimised Parabolic (TOP) nozzles, 
and have found use on a variety of actual launch vehicles because they perform better when 
over-expanded at ground-level altitude than the actual optimised bell nozzle (flow separation 
from the TOP nozzle wall is delayed at high back-pressure). 

Rao’s parabolas are known in Europe as quadratic Bezier curves, after French car bodywork 
draughtsman Pierre Bezier. Bezier curves are now used extensively in computer graphics. 

The shape of the bell nozzle changes only minutely with the propellants used (varying ratio of 
specific heats y) so one TOP nozzle methodology fits all propellants and is described below. 
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The geometric approach 


Construction 

The ‘Rao nozzle’ starts with Rao’s preferred throat geometry, presumably a Rockedyne 
(where he worked) optimisation. 



^ throat 


hozzle 

AXIS 


Two circular arcs are drawn: 


The first curve, of radius 1.5 R t , is drawn from an angle of say, -135 degrees, to the throat at 
-90 degrees (angles measured from the arc’s origin). Then the second curve of radius 0.382 
Rt is drawn from this angle of -90 degrees to an angle of (9 n - 90) at inflection point N. 

(Rt is the throat radius, R e is the exit radius.) 

Then a skewed parabola is drawn from point N to nozzle exit point E, tangent to the throat 
curve, and starting at an angle of 9 n and ending at an angle of 9 e . 


The radius of the nozzle exit: R e = \feR t equ. 2 


and nozzle length L N = 


V tan(15) ) 


equ. 3 for an 80% bell from equ. 1 


Angles 9 n and 9 e were pre-calculated by Roe to match his bell nozzle, and presented as 
graphical data from which the following chart is reproduced for various percent lengths (from 
Ref. 1): 
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(The data for expansion ratios greater than 50 is extrapolated). 

A parabola is then sketched-out using an ancient geometrical method for drawing a parabola 
detailed in Ref. 2. 



Straight lines are drawn at angles 9 n from point N, and 6 e back from point E, terminating 
where these lines cross at point Q. 

Next, both of these lines are divided into an equal number of divisions (four in this example) 
These are labelled a,b,c and e,f,g. 
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A straight line is then drawn from point a to point e, then from b to f, and c to g. These form a 
mesh, the edge of which gives the parabola outline. The parabola is also tangent to the lines 
QN and QE. 

Using many more divisions, on a CAD package for example, gives a sharper contour. 
Removing most of the mesh gives a series of straight-line segments: joining the midpoint of 
each segment with a smooth curve such as a CAD spline gives the nozzle contour. 



Alternately, interpolating along the construction lines gives the same points. In our example 
with four divisions and (4-1) construction lines, the parabola is defined at 1/4 the distance 
along line a-e, 2/4 along line b-f, and 3/4 along c-f. 

Efficiency 

At a length ratio of 85% bell, a nozzle efficiency of 99% is reached, and only 0.2% of 
additional performance can be gained by increasing the length ratio to 100%. For this reason, 
85% is often taken as upper bound. At length ratios below 70%, nozzle efficiency suffers. For 
these reasons, the 80% bell parabola is often chosen. 

Example 

An 80% bell nozzle with an area ratio of 70 (a typical upper stage nozzle) has angles 6 n = 33 s 
and 9 e = 7- from the above chart. Drawn on a CAD package it looks like this (mesh removed 
for clarity): 
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The mathematical approach 

The throat 

The equations of the above circular arcs defining the throat are defined trigonometrically, 
defining the origin of the coordinates as the centre of the narrowest part of the throat: 


For the entrant section: 
x = 1.5 R t cos 9 

y = 1.5 R t sin 6 + 1.5 R t + R t equ.s 4 

where: -135 < 9 < -90 

(The initial angle isn’t defined and is up to the 
combustion chamber designer, -135 degrees is 
typical.) 

For the exit section: 
x = 0.382 R t cos 9 

y = 0.382 R t sin 9 + 0.382 R t + R t equ.S 5 
where: -90 <9<(9 n - 90) 



The bell 

The bell is a quadratic Bezier curve, which has equations (see Wikipedia): 

x(t) = (1 — t) 2 N x + 2(1 — t)t Q x + t 2 E x 0 < t < 1 

y(t) = (1 - t) 2 N y + 2(1 - t)t Q y + t 2 E y 0 < t < 1 equ.s 6 

Selecting equally spaced divisions between 0 and 1 produces the points described earlier in 
the graphical method, for example 0.25, 0.5, and 0.75. 

Equations 6 are defined by points N, Q, and E (see the graphical method earlier for the 
locations of these points). 

Point N is defined by equations 5 setting the angle to (On- 90). 

Coordinate E x is defined by equation 3, and coordinate E y is defined by equation 2. 

Point Q is the intersection of the lines: Wq = m x x + C x and: QF = m 2 x + C 2 equ.s 7 
where: gradient m 1 = tan(0 n ) , gradient m 2 = tan (9 e ) equ.s 8 
and: intercept C 1 = N y - , intercept C 2 = E y - m 2 E x equ.s 9 
The intersection of these two lines (at point Q) is given by: 


Qx 


(c 2 -ci) n 

’ y 


(m 1 C2-m 2 Ci) 

(rrii-mz) 


equ.s 10 
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